Randomized Spectral and Fourier-Wavelet Methods for multidimensional Gaussian random vector fields

We develop some new variations of the Randomized Spectral Method (RSM) and Fourier-Wavelet Method (FWM) for the simulation of homogeneous Gaussian random vector fields based on plane wave decompositions. Some general conditions on the spectral density tensor of the random field which guarantee the theoretical convergence of the numerical representations are given. We then compare the various random field simulation methods through the examination of their efficiency and the quality of the Eulerian and Lagrangian statistical characteristics of a three-dimensional isotropic incompressible random field as simulated by ensemble or by spatial averaging. We find that the RSM, which is easier to implement, can simulate random fields with accurate second order Eulerian and Lagrangian correlations more efficiently than the FWM when the statistics are collected through ensemble This work is supported partly by the RFBR Grant 09-01-00152-a. ∗Corresponding author: Phone +001 (518) 276-6896, Fax +001 (518) 276-4824 Email addresses: [email protected] (Orazgeldi Kurbanmuradov), [email protected] (Karl Sabelfeld), [email protected] (Peter R. Kramer) Preprint submitted to Journal of Computational Physics February 5, 2012 averages. When the correlation structure is inferred from spatial averages of a single realization of the random field (through appeal to an ergodic theorem), the FWM is more efficient in simulating an accurate Eulerian correlation function, and both the FWM and RSM are of comparable efficiency in simulating second order Lagrangian statistics accurately.